Slow-Roll Inflation Preceded by a Topological Defect Phase a la Chaplygin Gas 
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We present a simple toy model corresponding to a network of frustrated topological defects of 
domain walls or cosmic strings that exist previous to the standard slow-roll inflationary era of the 
universe. Such a network (i) can produce a slower inflationary era than that of the standard scenario 
if it corresponds to a network of frustrated domain walls or (ii) can induce a vanishing universal 
acceleration; i.e., the universe would expand at a constant speed, if it corresponds to a network 
of frustrated cosmic strings red. Those features are phenomenologically modeled by a Chaplygin 
gas that can interpolate between a network of frustrated topological defects and a de Sitter-like or 
a power-law inflationary era. We show that this scenario can alleviate the quadruple anomaly of 
the cosmic microwave background spectrum. Using the method of the Bogoliubov coefficients, we 
obtain the spectrum of the gravitational waves as would be measured today for the whole range of 
frequencies. We comment on the possible detection of this spectrum by the planned detectors like 
BBO and DECIGO. 

PACS numbers: 98. 80. Bp, 04.30.-w, 95.36.-|-x 



I. INTRODUCTION 

The inflationary scenario explains several cosmological 
problems such as the spatial flatness and homogeneity 
of the universe, and it predicts an almost scale-invariant 
power spectrum of the primordial density perturbation, 
which agrees with most of the WMAP observation on 
cosmic microwave background (CMB) anisotropy. On 
the other hand, there remains some discrepancies be- 
tween the standard cosmological model predictions and 
the WMAP observations on large angular scales. One of 
them is the suppression of the I = 2 quadruple mode of 
the CMB temperature power spectrum (c.f. for example 
Q). The rc have been different approaches to solve this 
problem. In Refs 0,13, it is considered that the slow- 
roll inflation was preceded by a fast-roll stage where the 
kinetic energy of the scalar field was dominant, which 
leads to a suppression of the quadrupole moment. Such 
a fast-roll stage can be reached by tuning the initial con- 
dition of the inflation, for example. Another way to reach 
an early stage of kinetic energy dominance was consid- 
ered in Refs. 0, through bounces or cyclic universes 
where just before the bounce or the new cycle the kinetic 
energy density of the field dominates over its potential. 
Other authors have considered different pre-infiationary 
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phases of the universe to suppress the primordialpower 
spectrum at large scales. For example in Refs. p, 0]; a 
radiation-dominated pre-inflationary era was considered 
and in , a pre-inflation matter era provided by primor- 
dial micro black hole remnants was introduced. 

The CMB quadruple anomaly affects the lowest modes 
of the primordial power spectrum; i.e. modes with 
comoving wave numbers of the order of 10~^Mpc~^. 
Therefore, these are the flrst modes to exit the hori- 
zon and the last ones to reenter. Consequently we ex- 
pect that these modes are heavily affected by the physics 
of the very early universe, possibly the physics prior to 
the slow-roll inflationary era. Based on this reasoning, 
we propose a new cosmological period just before the 
slow-roll inflationary era. This period corresponds to an 
era described by a network of topological defects, which 
we will assume to be frustrated domain walls or cosmic 
strings , for the sake of simplicity. We expect that 

the production of topological defects at those scales fol- 
lows the prediction of high energy physics. In addition, 
given that the topological defects era proceeds the slow- 
roll inflationary era, the topological networks will affect 
mainly the lower modes of the power spectrum of the 
scalar and tcnsorial perturbations. Afterwards, they will 
soon be diluted during most of the inflationary era. 

There is a simple way to model the ideas previously 
mentioned. Wc consider the matter content of the early 
universe to be described by a kind of generalized Chaply- 
gin gas (GCG) [T^ - [l8| . which offers a smooth transition 
between different eras: 

• A period dominated by a network of frustrated do- 
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main walls (NFDW) at earlier time and a de Sitter- 
like phase at later time. 

• A period dominated by a network of frustrated cos- 
mic strings (NFCS) at earlier time and a de Sitter- 
like phase at later time. 

The idea of an early Chaplygingas phase in the universe 
was first suggested in Ref. [l9i (see also [l^, HH) and 
later extended in Refs. (23 - |24| . The energy density of a 
network of frustrated topological defects (NFTD) can be 
described in a compact way, for example, as 



P = 



B 



7,/3l(l+ai) 



-4i 



1/(1 + Ql) 



(1.1) 



where a is the scale factor, Bi and Ai are constants re- 
lated to the energy scale of the NFTD and the de Sitter- 
like inflationary era, respectively, ai and /3i are constants 
such that /?! = 1,2 for the network of domain walls and 
cosmic strings, respectively. We assume that 1 -|- ai is 
positive such that the inflationary era is preceded by 
a topological dominance phase. Let us be reminded in 
this regard that the energy density of NFDW and NFCS 
scales as 1/a and 1/a^, respectively |25|. It is worthy 
to stress that the NFTD epoch preceding the slow-roll 
inflationary can in principle produce inflation as well; 
indeed this is the case for NFDW, but this inflation is 
much slower, i.e. much lazier than the slow-roll inflation. 
Moreover, for a NFCS dominated period the universe is 
increasing its size at a constant speed; i.e. with no ac- 
celeration or deceleration. From now on whenever we 
refer to a pre-inflationary era, we will be referring to a 
pre-slow-roll inflationary era. 

This paper is organized as follows: In Sec. [ill in oi- 
der to solve the CMB quadrupole problem, we propose a 
NFTD era just before the standard inflationary era. We 
model this idea using the framework of matter a la Chap- 
lygin gas during the early universe as a way of interpolat- 
ing between a NFTD and a slow-roll inflationary period. 
We then explain how to constrain our model observation- 
ally. In Sec. mil we briefly review the standard cosmolog- 
ical perturbation formalism and obtain numerically the 
curvature power spectrum as well as the CMB tempera- 
ture anisotropy spectrum for the model presented in Sec. 
HTl In Sec. IIV| we obtain numerically the spectrum of the 
gravitational waves for the model introduced in Sec^Jlll 
by using the method of the Bogoliubov coefficient \2t 
Finally we summarize and conclude in Sec. |Vl 



II. MODEL BUILDING AND PARAMETERS 
FIXING 

We consider a spatially flat Friedmann-Lemaitre- 
Robertson- Walker (FLRW) universe filled with the mat- 
ter content described by Eq. (|l.ip . The energy conserva- 
tion gives 



where a dot corresponds to a derivative with respect to 
the cosmic time and H stands for the Hubble rate. By 
inserting Eq. (|l.ip into Eq. (|2.ip . one obtains the pressure 
of the matter content 
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Pi 



(2.2) 



In the Planck unit [c = h = G = \), the Friedmann 
equation reads 



(2.3) 



where = 8ttG = Stt. The conformal time r can be 
expressed as 

r = ^W"-' (^) " y^F.ic. l-b,c+l,y), (2.4) 

K C V^l / 

where y = [l + (Bi/Ai) a->^^^^+°'^^~\ b = 
l/[/?i(l + ai)], c = (-2 + /3i)Ml + ai)/3i], and 
2F1 is a hypergeometric function [23|. From Eq. (|2.4p 
we can see that the universe began in NFTD dominated 
era in the past infinity, and turned into a de Sitter-like 
space, in which a oc 1/t, at later time. 

In order to obtain the power spectrum of the gravi- 
tational wave as would be measured today, which will 
be obtained in section IIVI we divide the expansion of 
universe into three successive periods: the pre-inflation 
NFTD dominated era, the inflating phase, and the stan- 
dard ACDM epoch. The energy density of each of these 
periods can be modeled as 
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(2.5a) 

(2.5b) 
(2.5c) 



p + 3i/(p + p) =0, 



(2.1) 



Expression (|2.5aP describes the energy density of the 
NFTD era, which was introduced in Eq. (jl.ip . followed by 
the de Sitter-like inflating phase. The model described 
by Eq. (|2.5b[) was previously studied within an infla- 
tionary framework in Ref. [22j (see also Ref. [l^) and, 
under suitable constraints on A2, B2, and a2 , can depict 
the transition from the de Sitter-like era to the radiation 
dominated era. The energy density (|2.5cp is the standard 
ACDM model, in which piOi PmOi and p\ are the energy 
densities of the radiation, matter, and dark energy to- 
day, respectively. As we will show later, the parameters 
of the model can be constrained using observational data 
corresponding to the present energy density of radiation, 
the scalar power spectrum, and the spectral index at a 
given pivot scale. In addition, by requiring that the en- 
ergy density is continuous at each transition, we have the 
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conditions 



So 



4\ 1 + Ct2 



(2.6) 
(2.7) 



In order to obtain the scalar power spectrum, it is use- 
ful to model the matter content in the first two periods of 
Eq. ((231) : i.e. those described by Eqs. ((23a)) and (j2.5bp . 
through a scalar field, with the condition that the energy 
density and pressure of the scalar field are the same as 
that given by Eqs. (|2.5ap and (j2.5b[) . The energy density 
and pressure of the scalar field are 



P0 



(2, 



where the primes denote the derivatives with respect to 
the conformal time. The potential of the scalar field in 
the first period is given by 



^i(0) = Y](6-/3i)cosh 



K (1 + ai) Vpi 2 



l + ai 



+ Pi cosh 



K (1 + ai) VPi-^ 



l + Ql 



(2.9) 



where Vq 



ll/(l+ai) 



Similarly, the potential in the 



second period can be obtained by replacing ai with 02 
and setting /3i = 4 in Eq. (|2.9p . which gives 



1 + Q2 



1^2(</') = y {cosh[K(l + a2)< 

+2 cosh [k (1 + 02) 0] ^ 



(2.10) 



The result coincides with the potential in Ref. [22|, as it 
should be. 

Next we tackle the issue of analyzing the potentials 
(|2.9p and (|2.10p . First of all, we consider that the scalar 
field potential p.9p has a unique minimum at (/> = to 
maximize the amount of infiation during the first period. 
Notice that unless this condition is imposed, the scalar 
field might roll down the potential till it reaches the min- 
imum of Vi{(f>) and then would have to climb up to reach 
the local maximum located at = 0, as shown in Fig. [1] 
Imposing that the potential (j2.9p has a unique minimum 
reached at = implies a condition on the parameters 
ai and /3i that 



ai < 



6 -/?i 

(3i 



(2.11) 



Therefore, bearing in mind that (i) (3i — 1,2 for NFDW 
and NFCS, respectively, and (ii) < 1 4- ai so that the 
phase of FNTD precedes the inflationary phase, we con- 
clude that -1 < ai < 5 for NFDW and -1 < ai < 2 
for NFCS. We show the shape of the potential Vi{(j)) for 
different cases when the condition (|2.1ip is fulfilled and 
violated in Fig. [TJ 
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FIG. 1. The scalar field potential for the NFDW era for dif- 
ferent values of ai. The solid curve corresponds to qi = 1 and 
the dashed curve corresponds to ai—7. It is similar for the 
case of NFCS, which we omit for simplicity. 



In addition, we can constrain our model; i.e. the poten- 
tials p.9p an d ([2.10p . using the methodology used by one 
of us in Ref. [22] . More precisely, we can use the WMAP7 
observation of the power spectrum of the comoving cur- 
vature perturbation, Pg = 2.45 x 10~^, and the spectral 
index, Ug = 0.963, at the pivot scale fco = 0.002 Mpc~^ 
to fix the parameters in our model [1^. We can as well 
impose a bound on the number of e- folds, Nc, since a 
given mode exits the horizon until the end of inflation 
as done in Ref. [22[. This gives the best- fit values for 
a2, Vq, and, therefore, A2. Notice that once Vq is fixed, 
the parameter Ai is fixed for a given ai as well, since 

The parameter Bi in Eq. (|2.5aP fixes the energy density 
of the NFTD, which strongly affects the lowest modes 
that exited the horizon around the onset of inflation, 
and causes a significant drop on the lowest modes of the 
primordial spectrum of the curvature perturbation. Al- 
though we expect that the NFTD would affect the lowest 
modes, we must make sure that the observed curvature 
power spectrum Pg and the spectral index Us at the pivot 
scale ko are recovered. Therefore, we choose the value of 
Bi such that Pg and Ug correspond to the correct values 
at ko, and the amplitude of Pg drops at the scales whose 
comoving wave numbers are smaller than kg. Roughly 
speaking, the parameter Bi controls the horizontal shift 
of the curvature power spectrum. 

However, following this procedure, it turns out that 
we can not find a set of values for the parameter Bi that 
satisfies the constraint on the spectral index. In fact, 
in this model Bi turns out to be always smaller than 
0.9. The reason of this drawback is the second period 
described by Eq. (j2.5bp . which corresponds to a short 
transition from a de Sitter era to a radiation dominated 
period. More precisely, the model described by Eq. (|2.5p 
does not give enough e-folds during the inflationary era. 
We show, as an example, in Fig. [2] how the scalar field 
rolls too quickly and the radiation dominated phase is 
reached too early in the case corresponding to a NFDW. 

We therefore suggest an alternative model described 
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FIG. 2. The scalar field potential for /3i = 1. The upward 
concave curve is Vi (cf. Eq. (|2.9p ) and the downward convex 
curve is V2 (cf. Eq. (I2.10|l '). We can choose (j) < ior Vi and 
> for V2 to describe the potential of the scalar field 0. 
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A2 



A. 



, B2 



ji+ft a4(i+/33) 



(2.12a) 



(2.12b) 



fto(^)^ + Pmo(^)-VpA, (2.12c) 

where /3i = 1, 2 discriminates the NFDW and the NFCS 
as in the former model, and /32, P3, Bi, A2, B2 are con- 
stants we will explain below. We choose this model for 
the second period since it generates an almost flat cur- 
vature spectrum for modes larger than the pivot scale 
ko = 0.002 Mpc~^ and gives enough e-folds during the 
power law inflationary period. In addition, the model 
introduces in a natural way that a NFTD precedes the 
power law inflation as (1 + /32)/(l + /^s) < /3i (please see 
also the conditions ((2?T4l) . 0J5\i and ((2T5|l '). 

The first period in Eq. p.l2ap describes the matter 
content of the universe during a period that transits from 
a NFTD dominated phase to a power-law inflationary 
era. The parameters Bi and A2 are associated with the 
energy scale of the NFTD and the power-law inflation, 
respectively. The second period with the energy density 
(j2.12b[) was previously studied within another inflation- 
ary framework in Ref. (23j . It connects smoothly a power- 
law inflating phase with a radiation dominated universe, 
and the constraints on the parameters and arej^ 



1 + /32 < 0, 
1 + /33 < 0, 
2(l + /33) < 1 + /32. 



(2.13) 
(2.14) 
(2.15) 



^ The notation is different from the one used in the work l23l . The 
parameters /3 and a in Ref. |23|| are denoted as /32 and here, 
respectively. 



These constraints imply that (i) there is a power-law in- 
flating phase, (ii) the inflationary era precedes the radia- 
tion dominated period, and (iii) the null energy condition 
is always fulfilled so that there is no super-inflationary 
phase. Finally, the energy density described by (|2.12c[) 
corresponds to the ACDM model as that described by 

Eq. dm . 

Although there seems to be many free parameters, they 
can be fixed down to only one by the following procedure: 
(i) Fix B2 by the current amount of radiation for a given 
(3^. (ii) Constrain the power-law expansion quantified 
by ^(i+/33)/(i+^'2) ^j^^ WMAP7 data of the curvature 
power spectrum Ps and the spectral index n^. Since there 
arc three parameters {A2, P2 and P^) to be constrained 
by only two conditions (Pg and n^), we are left with the 
only one free parameter which we choose to be /33 for 
convenience, (iii) Fix Bi such that Pg and rig give the 
correct value at the pivot scale fco and the drop on Pg 
just begins at comoving wave numbers smaller than fco. 

Again, it is suitable to introduce a scalar field that 
mimics the matter content described in Eqs. (|2.12ap and 
(|2.12bP : i.e. we describe the dynamics of the model 
through a scalar field with a potential whose energy den- 
sity and pressure can be obtained from Eq. p.Sp . During 
the NFTD period (cf. Eq. (j2.12ap ). the mapping between 
the scalar field and the perfect fluid of our model leads 
to 



0(a) = \/v coth ^ \ 



/Atanh-\/c+(l-C)^, (2.16) 



Vi{a) 



(6-/3] 



A. 



(2.17) 



where C = [l + {Al,''''^^" / Bi)a''^--]-\ v = {l+p2)/il + 
Ps), and Vi(0) stands for the scalar field potential during 
this period. Similarly, we map the perfect fluid with the 
energy density (cf. Eq. (|2.12bp ) to the continuation of 
the scalar field (p with a new potential V2 {4>) [1^ 



0(a) = - 
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FIG. 3. This plot shows the rescaled potentials given in 
Eqs. (|2.16|) - (|2.19[) versus the scalar field 4>- The blue curve 
corresponds to the NFCS case. The red curve corresponds to 

the NFDW case and Vb = Ay^^+^^^ (^2/52)"^'+''=^'''^''^'+''='". 
The energy scale of inflation, Vb, in our model is about 10^^ 
GeV for both NFDW and NFCS. 



where ^ = (53/^2)0'^ and q = l + /32-4(l + /33). The po- 
tential (|2.19p was previously obtained in Ref. [IJ. Such 
a potential, with an appropriate initial condition, drives 
a power law inflation and mimics a radiation dominated 
universe afterwards. 

It is worthy to notice that an integration constant ap- 
pears when we integrate Eq. (j2.8p after mapping it to 
the energy density and pressure of a given perfect fluid. 
Therefore, we can always choose the constant properly 
such that the potentials Vi((/)) and V2{4>) are continuous 
at the connecting point. We show an example of the po- 
tentials Vi[(j)) and V2((/)) in Fig.[3l The scalar field starts 
with a negative value and rolls down the potential as the 
universe inflates until it reaches the radiation dominated 
era. 



III. SCALAR PERTURBATIONS 

The accelerated expansion of the universe during the 
primordial inflationary era converts the initial quantum 
fluctuations in the universe into macroscopic cosmologi- 
cal perturbations, which leads to the inhomogeneity we 
observe nowadays in the CMB [2^ [sO]. Following the 
standard approach, we will use gauge invariant quanti- 
ties that involve the metric perturbations and the scalar 
fleld fluctuations [3l[. For convenience, we will choose 
the comoving curvature perturbation, R, which in addi- 
tion is conserved on large scales [H, [s^ . 

We expect that a NFTD in the very early universe and 
just before the inflationary era could give the appropriate 
corrections to the quadrupole modes of the CMB data as 
observed nowadays [28j . We will next quantify the quan- 
tum cosmological perturbations during that period and 
obtain the power spectrum of the scalar perturbations. 



A. Formalism 

The scalar perturbations can be described by introduc- 
ing the variable (see, for example, Ref. [sij ) 



(3.1) 



where z = The variable u can be decomposed into 
Fourier modes, Uk, which fulflll the field equation [s^ ] 



1 (fz 
z dr^ 



Uk = 0. 



(3.2) 



The modes Uk can be mapped to the spectrum of the 
comoving curvature perturbations which reads [13] 



27r2 

/c3 



■PR{k) 



\Uk\ 



(3.3) 



Given that we are dealing with adiabatic perturbations, 
the comoving curvature perturbations remain constant 
on large scales and consequently we can equate the power 
spectrum at the horizon exit with the power spectrum of 
the primordial scalar perturbations at the horizon reentry 
as observed on the CMB. Therefore, for a given mode fc, 
the spectrum is evaluated at the horizon exit; i.e. when 
k = QcmssH , where across stands for the value of the scale 
factor when the mode exists the horizon. 



B. Numerical Solutions 

We next obtain the evolution of the mode function 
Uk (r) for each comoving wave number k in order to obtain 
the curvature perturbation spectrum. We will tackle this 
issue numerically rather than using the standard results 
for slow- roll inflation [2^ [13, [HI, because those condi- 
tions are not fulfilled at very early time when the NFTD 
is dominant. It is easier to solve Eq. (|3.2p numerically by 
splitting it into two first order differential equations, 




(3.4) 



where we have set X = Uk- 

In addition, we need to impose a set of boundary condi- 
tions at the time when the wavelength of a given mode k 
is much smaller than the Hubble radius; that is, k ^ aH. 
Those boundary conditions will depend on the specific 
NFTD scenario and also on the given scale of the mode, 
as wc will explain shortly. It is also worthy to stress two 
things: (i) the approximation z" / z « a" /a holds during 
the NFTD dominance and during the power law infla- 
tionary era (See Fig. |4] and Fig. [5]), and (ii) the analytical 
solutions for the power-law inflation [36| can be used as 
boundary conditions for the modes with comoving wave 
numbers larger than roughly 10"'^ Mpc~^ as explained 
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FIG. 4. The green dashed curve corresponds to z" j z and the 
black curve corresponds to a" /a. As can be seen that the ap- 
proximation z" / z ~ a" /a holds during the NFDW dominance 
and the power law inflationary era. 
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FIG. 5. The green dashed curve corresponds to z" / z and the 
black curve corresponds to a" /a. As can be noticed that the 
approximation z" / z ~ a" /a holds during the NFCS domi- 
nance and the power law inflationary era. 



next. We first remind that the power law solutions for X 
and Y read ^ 



(3.5) 



(3.6) 

where I is the exponent characterizing the power law ex- 
pansion in terms of the conformal time r; i.e. a (x tK 
For our model, / = ((1 + /^2)/(2(l + Ps)) For 
those modes whose k > lO^'^Mpc"^, we can start the 
numerical integration of Eq. (|3.4p during the power-law 
inflationary dominated era where the condition k ^ aH 
is still satisfied. However, for scales roughly smaller than 
10"'^ Mpc~^, we split the boundary conditions imposed 
on the NFDW and the NFCS separately. Despite the 
general solutions (|3.5p and (|3.6p are still valid for the 
NFTD, the exponent describing the power law expan- 
sion, I, depends on the specific characters of the NFTD: 



• During NFDW dominant era, p ^ Bi/a^^ with 
/3i = 1, which also implies a power-law expansion 
because a(r) cx T^/'^fi^/'^~'^) . Thus we can use the 
solutions p.Sp and p.6p as boundary conditions 
with I = l/(^i/2- 1). 

• During the NFCS dominant era, p ^ B\ja^^ with 
/3i = 2. Notice that for this value of /3i the expo- 
nent / in Eqs. p.5p and p.6p is not well-defined, so 
we need to fix the initial condition in this case in a 
different way. It can be shown from Eq. (j2.3p that 
a"/a is a constant if p oc Therefore, we can 
define a new variable fc' = y^/c^ — a" j a and inter- 
pret k' as an effective comoving wave number. The 
wave equation (|3.2p then becomes 



dr2 



(3.7) 



which has two properly normalized linearly inde- 
pendent solutions. 



Mfc(r, k) 



(3.8) 



We choose the exponent with the minus sign be- 
cause it reduces to the Minkowski initial condition 

m. 

Finally, it is easier to use the scale factor as the inde- 
pendent variable in the numerical integration of Eq. (j3.4p 
instead of the conformal time. The relation between the 
conformal time t and the scale factor a is 



V3 



1 



:Bi 



-1/2 



K /l(/?l-C) 

2F1 [/i,l-.g,/i + l; 



Bi 



(3.9) 



-133 



'/B,)a 



/3i-Cl-l 



g 



(1 



where a; = 1 — [1 

l/(2C))/(/3i - C), and e = (/3i/2 - l)/(/3i - C)- 

The resulting curvature power spectra are shown in 
Fig. [HI and Fig. [71 corresponding to the NFDW and the 
NFCS, respectively. Let us recall that all the parame- 
ters used here have been fixed by imposing the obser- 
vational constraints in the way stated in Sec. |TT1 For 
modes k such that 10~^Mpc-^ < k < lO^Mpc-^, we 
obtain a constant slope on the logarithmic amplitude of 
the curvature perturbation. This is a simple consequence 
of the intermediate phase given in Eq. (j2.12bp previously 
analyzed in Ref. [23| , and implies a power- law inflation. 
Our new and important result is the drop of Pr for the 
modes whose k <10~^Mpc~^, which is helpful in explain- 
ing the quadrupolc anomaly through an alternative way 
from those used in Refs. 0, S S Ifl Such a decrease 
of Pr is a consequence of the NFTD era just before the 
inflation. 

We calculate the CMB temperature anisotropy spec- 
trum by the numerical package CMBFAST j37H39j with 
minor modiflcations to the form of specifying the primor- 
dial power spectrum. In the modifled version of CMB- 
FAST, the primordial power spectrum is fed into the code 



log[P«] 
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FIG. 6. This plot corresponds to the curvature perturbation 
spectrum for the inspired modified GCG model with /3i = f 
(see Eq. (I2.f2a|l 'l. which describes a NFDW dominated era 
followed by a power-law inflationary period. We choose ~ 
— f .05. The vertical dashed line corresponds to the pivot scale 
k = 0.002 Mpc-\ 




Muliipok Mumeni (1) 

FIG. 8. The CMB temperature anisotropy spectrum. The 
dots with error bars are the WMAP 7-year data. The solid 
line is the prediction of standard inflation with a power law 
spectrum. The dashed and dotted lines are the spectra of the 
scenarios of NFCS and NFDW, respectively. 




FIG. 7. This plot corresponds to the curvature perturba- 
tion spectrum for the modifled GCG model with jSi — 2 
(see Eq. (|2.12ap '). which describes a NFCS dominated era 
followed by a power-law inflationary period. We choose 
Pa = —1.05. The vertical dashed line corresponds to the 
pivot scale k = 0.002 Mpc"^ 



IV. TENSOR PERTURBATION 

As it is well known within the inflationary paradigm, 
aside from the scalar perturbations, there are also ten- 
sorial perturbations which give rise to a stochastic back- 
ground of gravitational waves generated from the vacuum 
quantum fluctuations during the inflationary era. If the 
stochastic background of GWs of cosmological origin is 
ever detected (hopefully on the next decades [10, |4l| ) , it 
will provide an amazing amount of information and ad- 
ditional probes about the very early universe. We next 
obtain the GWs power spectra predicted by the model 
described in Eq. (|2.12|) . Our analysis will be based on 
the method of the Bogoliubov coefficients and this study 
is motivated by the fact that we hope the detection of 
the GWs power spectrum can throw some light on the 
physics of the early universe, as we have already men- 
tioned. 



as an interpolating function instead of a functional form. 
This adjustment is made so that rather than accept- 
ing only the nearly-scale-invariant spectrum as its initial 
condition, CMBFAST is now compatible with any gen- 
eral shape of initial spectrum. This feature is essential 
to our case since the primordial spectra obtained from 
our scenarios severely deviate from the scale-invariant 
form in large scales whose wave numbers are smaller than 
10-3Mpc"\ 

Fig. [8] shows the CMB spectra generated by the NFTD 
scenarios along with the standard inflationary model. 
The data of WMAP 7-year observation [2^ are also 
marked in the plot. Note that the NFCS has a stronger 
effect on reducing the amplitude of the lower modes than 
the NFDW does. This is a result of the initial slope and 
the turn-around point of the curvature power spectrum 
induced by NFDW and NFCS as shown in Fig. |1] and [T] 



Formalism 



The gravitational waves energy density is given by [42| 



PGW ^ J (Nkir)) {2f^)^dLU, (4.1) 

where (Affc(r)) is the number of gravitons created for a 
given mode k in a different vacuum state as the universe 
expands. Then the dimensionless relative logarithmic en- 
erg y sp ectrum ficw a-t the present time tq can be defined 
as [42] 

S2gw(w,To) = — — r-7j (to) 



Pc(to) dlnw 



(4.2) 
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We will use the method of the Bogoliubov coefficients to 
obtain the number of gravitons created as the universe 
expands [26| . 

The tensor perturbation hij can be decomposed into 
Fourier modes which obey: 



dr2 



/2 Id^aV ^ 



(4.3) 



The initial vacuum state changes as the universe ex- 
pands. This change can be described through a Bogoli- 
ubov transformation: 



aik)^akA{k)+l3kA\~k), 



(4.4) 



where ak and f3k are the Bogoliubov coefficients. In par- 
ticular, the coefficient /3fc gives the number of gravitons 
created for each mode k; i.e. (A^fe(r)) = |/3fe(r)|. Given 
that the universe was initially in a vacuum state, we im- 
pose the initial conditions that = 1 and /3fc = 0, which 
corresponds to the absence of gravitons at the beginning. 
In addition, A and A^ correspond to the annihilation and 
creation operators of the gravitons. 



B. Numerical Solutions 

The evolution equation for the Bogoliubov coefficients 
can be simplified through the substitutions [2^ 



X{k;T)+Y{k;T) 

26 (r) 
X(fc;T)-y(fc;T) 



(4.5) 
(4.6) 



where ^ is the appropriate boundary condition chosen 
for a given mode k. A discussion parallel to that given 
in Sec. IIIII follows: 

• For the modes with comoving wave numbers 
roughly smaller than 10^'^ Mpc^^, the integration 
starts during the NFTD dominated era where the 
condition k 3> aH is satisfied and where the mode 
is well inside the horizon. Therefore, we can just 
use the Minkowski initial condition 12911 




10-1° 100 10'"^ 10^^ 10^^ 10^" 
a 

FIG. 9. The integration for a given mode k. The dashed curve 
corresponds to (aH)^ and the solid curve corresponds to a" /a 
for a modified /3i = 1 GCG model. The horizontal red Une 
shown in the figure corresponds to the pivot scale k = 0.002 
Mpc~^, and the horizontal dashed blue line corresponds to 
the intersection of the first two periods. 




38 



,50 



lO"'*' 100 10^4 10^^ lO-^'* 10 
a 

FIG. 10. The integration for a given mode k. The dashed 
curve corresponds to (aH)^ and the solid curve corresponds 
to a" /a for a modified Pi — 2 GCG model. The horizontal 
red line shown in the figure corresponds to the pivot scale k = 
0.002 Mpc"'^, and the horizontal dashed blue line corresponds 
to the intersection of the first two periods. 



(4.7) 



for the NFDW period. For the NFCS, we use the 
initial condition given in Eq. p.Sp . 



• For the modes with comoving wave number roughly 
larger than 10""^ Mpc~^, we can start the integra- 
tion during the power-law inflationary dominated 
era where the condition k ^ aH is satisfied and 
where the mode is well inside the horizon. Thus 
we can also use the Minkowski initial condition in 
Eq. (ITTl) . 



By substituting Eqs. ([45|) - (|4^ into Eq. ([Ol) and us- 
ing the initial conditions of ^ chosen above, Eq. (|4.3|) 
can be simplified to a set of first order linear differential 
equations of the variables X and Y [2^ 




0. 



(4.8) 



After integrating Eq. (|4.8p . we use Eq. (j4.6p to obtain 
the number of gravitons created till the present time tq 



{Nkir)) = \l3k(To)\ = 



\X-Y\ 



1 



(4.9) 
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FIG. 12. This plot corresponds to the gravitational wave spec- 
FIG. 11. This plot corresponds to the gravitational wave spec- tru m for th e inspired modified GCG model with /3i = 2 (see 
trum for th e inspired modified GCG model with /3i = 1 (see gq. pT2i)l l. which describes a NFCS dominated era foUowed 
Eq. (|2.12a[)), which describes a NFDW dominated era fol- by a power-law infiationary period with fis = -1.05. 
lowed by a power-law inflationary period with /Ss = —1.05. 



SUMMARY 



The integration method can be read from Fig. [5] and 
Fig. [To] for NFDW and NFCS, respectively. For a given 
mode fc, we start the integration at the time when the 
scale factor a is a thousandth of the scale where the mode 
exits the horizon; i.e. Ostart = O.OOlflcrossi- The integra- 
tion ends at the time when the scale factor a is thousand 
times of the scale where the mode reenters the horizon; 
i.e. acnd = lOOOacrossa- The reason for stopping our nu- 
merical calculations after the modes are well inside the 
horizon and not continuing the numerics till the present 
time is the following: once a mode is inside the horizon, 
the mode begins to oscillate and the number of gravi- 
tons created will not change anymore. Therefore, we can 
save a lot of computational time with this methodology 
[2^ . The resulting spectra are shown in Fig. [TT] and Fig. 
1121 The frequencies range approximately from 10^^^ to 
10'' Hz, and the dimensionless relative logarithmic energy 
spectrum of the GWs falls approximately between 10^^'' 
and 10""'^^. This result fulfills the current observational 
upper limits (all solid curves for various gravitational- 
wave detectors) shown in Fig. 2 of Ref. [i^. The gravita- 
tional waves predicted in our model could be detected in 
future experiments by some planned detectors like BBO 
and DECIGO whose sensitivities are shown in Fig. 6 of 
Ref. HH. It was shown in Ref. [H (see also Ref. [HI) 



that only large frequencies are sensitive to the slow-roll 
regime controlled by the parameter P^. For those fre- 
quencies, a larger gives rise to a higher maximum of 
the potential a" /a, a fact that translates into a shift of 
the power spectrum of the GWs towards the right. For 
middle frequencies, the plateau of the spectrum merges 
for different /S^ since the energy scale of inflation is al- 
most fixed and independent of the values acquired by /Ss. 
Similarly, we have shown that for low frequencies, the 
spectrum merges since the energy scale of the NFTD is 
almost indifferent to the changes of the parameter /S^ . 



The goal of this paper is to propose a possible solution 
to the low quadruple of the CMB which is based on a 
NFTD period just before the slow-roll inflationary era. 
The NFTD can be of two kinds: (i) a NFDW implying a 
slower inflationary rate than that implied by the standard 
slow-roll inflation, or (ii) a NFCS describing a universe 
sitting on the border line between acceleration and de- 
celeration. We model such a matter content for the very 
early universe as a kind of a generalized Chaplygin gas 
described by Eq. (|2.12p (see also Eq. (|l.ip ). which gives 
a smooth transition from a NFTD dominated era to a de 
Sitter-like phase or a power-law inflationary era. We con- 
strain our model using the WMAP7 data for the power 
spectrum of the scalar perturbations, P, = 2.45 x 10~^, 
and the spectral index, rig = 0.963, at the pivot scale, 
fco = 0.002 Mpc~^ HI]. After fixing the parameters of 
the model and imposing suitable boundary conditions for 
the perturbations, we obtain the curvature power spectra 
for the cases of a NFDW and a NFCS, as shown in Figs. [6] 
and[71 The most important feature of the spectrum is the 
drop of the power starting from fco ^ 0.002 Mpc~^. The 
reason behind this drop is the high energy effects of the 
NFTD on the lowest modes. Therefore, this mechanism 
can alleviate the observed low CMB quadrupole mode 
0, H, as shown in Fig. [8l Besides, we obtain the 

power spectra of the gravitational wave predicted by our 
model, which are consistent with the tensor-to-scalar ra- 
tio of the current observation ji^, |4l| . The regime of fre- 
quencies roughly larger than 10~^ Hz is within the reach 
of planned detectors such as BBO and DECIGO. Those 
frequencies are sensitive to the only degree of freedom of 
our model described by the parameter As a result, 
could be constrained by the projected sensitivities of 
BBO and DECIGO. 

In summary, even though the model presented here 
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could be seen as an over simplified toy model, it can give 
a guiding light towards a more consistent picture of a 
topological defect era prior to the slow-roll inflationary 
one, which could alleviate the observed CMB anomalies 
such as the suppression of the I = 2 quadruple mode of 
the CMB. 
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